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Abstract — In this paper, we propose a distributed Alamouti 
space-time code (DASTC) for two-way relay networks employing 
a single amplify-and-forward (AF) relay. We first derive closed- 
form expressions for the approximated average sum-rate of 
the proposed DASTC scheme. Our analysis is validated by 
a comparison against the results of Monte-Carlo simulations. 
Numerical results verify the effectiveness of our proposed scheme 
over the conventional DASTC with one-way communication. 

I. Introduction 

Relay networking has been considered as an efficient 
approach to increase the communication range of wireless 
systems. However, dual-hop half-duplex relay networks lose 
half of the throughput compared to the direct communication 
due to the fact that the relay cannot transmit and receive 
simultaneously. To overcome this drawback, a two-way (or 
bi-directional) relay network has been presented in Q], where 
two nodes, namely Si and S2, transmit simultaneously to the 
relay R in the first hop, and in the second hop the relay R 
forwards its received signals to both terminals Si and S2. With 
this strategy, this loss can be remarkably compensated. 

As a result, two-way relay networks have gained great 
attention in the research community (e.g., see Hl-flSI]). In view 
of the related work of two-way relay networks, the upper and 
lower bounds for average sum-rate have been investigated in 
|fl~). The performance of two-way amplify-and-forward (AF) 
relay with network coding has been investigated in [2|. By 
giving up the strict separation of downlink and uplink signals 
through either time or frequency division duplex, a two-hop 
relaying, namely space division duplex relaying, is proposed 
in 0. The exact closed-form expressions of error probability, 
average sum-rate for two-way AF relays have been presented 
in H, 0. 

Distributed Alamouti space-time code (DASTC) with AF 
relays originally applied for one-way relay networks ||6| has 
recently extended to two-way systems (see, e.g., 0-||9)). 
In particular, a two-way relaying scheme where two sources 
equipped with two antennas transmit Alamouti code through 
the help of an AF relay has been proposed in Q, 0. Upper 
and lower bounds of average sum-rate have also been derived 
for this particular bi-directional relaying system. Another 
system in which the burden of deploying multiple antennas 



on sources is transferred to the relay has been proposed in 
|9|. By assuming that the relay is equipped with two antennas 
and each source is equipped with a single antenna, an upper 
bound of symbol error probability has been obtained. As 
can be observed from all the above schemes, either relay or 
source requires multiple antennas which may be prohibited in 
practical implementations due to the high demand for low-cost 
and small-size portable devices. 

To alleviate this requirement and make the system realistic, 
a DASTC scheme which all terminals having only one antenna 
has been presented for one-way AF relay networks ifTOl . ATI . 
However, this system still faces loss in spectral efficiency. 
Hence, in this paper, we propose a DASTC scheme for 
two-way AF relay networks which significantly recovers this 
loss. More importantly, unlike the analysis work in 0, J8J 
where the bounds of average sum-rate have been shown, we 
have derived an asymptotically tight approximation for the 
average sum-rate of the proposed two-way scheme. The final 
expression is given in the form of Fox's H-function which 
enables us to investigate the performance of the proposed 
scheme. In addition, we also provide the numerical results 
to verify the correctness of our analysis. 

The rest of this paper is organized as follows. In Section fTTl 
we introduce the system model of the proposed DASTC for 
two-way AF relay networks. Then, in SectionUm we derive the 
tight approximation for the average sum-rate of the proposed 
scheme. Numerical results are shown in Section|lV]to validate 
the analysis. Finally, we conclude the paper in Section [V] 

Notation: Throughout the paper, we shall use the following 
notation. Vector is written as bold lower case letter and matrix 
is written as bold upper case letter. The superscripts * and 
f stand for the complex conjugate and transpose conjugate, 
respectively. I n represents the n x n identity matrix. \\A\\ F 
denotes Frobenius norm of the matrix A and \x\ indicates the 
envelope of x. E x {.} is the expectation operator over the ran- 
dom variable x. A complex Gaussian distribution with mean 
H and variance a 2 is denoted by C/V(/j,, c 2 ). Let us denote 
Af m (m, E) as a complex Gaussian random vector with mean 
vector m and covariance matrix S. T (a, x) is the incomplete 
gamma function defined as T(a,x) = t a ~ 1 e~ t dt and 
K, n (.) is the nth-order modified Bessel function of the second 
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Fig. 1. System model. 
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II. System Model of Distributed Alamouti 
Space-Time Codes with Two- Way AF Relaying 

In this section, we introduce the system model of the 
proposed DASTC scheme with two-way AF relaying. Consider 
a wireless AF relay network consisting of three terminals as 
shown in Fig. [1] Two sources, namely Si and S2, exchange 
the information through the assistance of an AF relay R. Each 
terminal is equipped with a single antenna which operates 
in half-duplex mode. We also assume that the channel for 
all links induces quasi-static fading, i.e., the channel remains 
constant for a block spanning over at least six symbols- 
intervals and varies independently for every block. Due to the 
channel reciprocity, we note that channel gains of A — > B 
and B — > A links, with A, Be {Si, S2, R}, are identical. 

The communications of the proposed two-way AF DASTC 
scheme occur in three phases. In the first phase, both Si and 
S2 transmit the first row of Alamouti code lfl2ll to R with the 
same transmit power per symbol V s . The received signal at 
the relay R is given by 



y R = hxSi + h 2 xi + n R 



(1) 



where s± = [ s\ s 2 ] and x\ = [ X\ x 2 ] are the two 
first rows of the Alamouti codes generated by Si and S2, 
respectively. The channel coefficients for the link from Si and 
S2 to R, respectively denoted as hi and h 2 , follow a Rayleigh 
distribution, i.e., hi — CA/"(0,Sli) and h 2 ~ CAf (0, fi 2 ). 
The vector n R is the complex additive white Gaussian noise 
(AWGN) with zero mean and variance No. 

In the second phase, R amplifies the received signal, i.e., 
y R , with the same power constraint as in the first phase and 
forwards to S2 while Si sends the second row of the Alamouti 
code to S2. The received signal at the source S2, y 2 , i s written 
as 



y 2 = h 2 Gy R + h s 2 + n 2 



(2) 



where G is the scalar amplifying gain at the relay, ho ~ 
CN (0, Hq) is the channel gain for the link between Si and S2, 
s 2 = [ — si } is the second row of the Alamouti code, 
and n 2 is the AWGN vector with zero mean and variance Nq. 



Similarly, in the third phase, R transmit Gy R whereas S2 
sends the second row of the Alamouti code to Si. The received 
signal at the source Si, y v is shown as 



J/i = hiGy R + h a x 2 +m 



(3) 



where x 2 



is the second row of Alamouti code 



and Hi is the AWGN vector with zero mean and variance N . 

The amplifying gain G is determined to satisfy the aver- 
age power constraint between the relay and source, i.e., R 
consumes the same amount of power as each source. In this 
paper, we assume channel state information (CSI)-assisted AF 
relay, i.e., the relay perfectly knows instantaneous values of 
the channel gains hi and h 2 . Taking these information into 
account, we have 



E 



{\\Gy R \\ 2 F }=E{\\ Sl \\l}=E{\\x t \\l} (4) 



yielding 



G 2 = \hi 



,2_ M, ,2 , , h |2 



7 



(5) 



where 7 = Vg/No is the average signal-to-noise ratio (SNR). 
With sufficient large SNR, the amplifying gain can be tightly 
approximated as 



(6) 



It is assumed that the two sources have perfect knowledge 
of the corresponding channel coefficients to fully cancel self- 
interference. Hence, the received signal at Si and S2 given in 
(f3]l and (|2]) can readily be formed, respectively, as follows: 



where h = 
Gh 2 n R + n 2 , 



yi = hX + zi 
y 2 =hS + z 2 

Ghih 2 ho ], Z\ = Ghin R 



(7) 
(8) 



S = 



si s 2 



and X 



- n u z 2 

X\ x 2 



Here, the new noise vector Zi, with i£ {1, 2}, is the complex 
Gaussian vector, i.e., z l — Kf 2 (o, N (l + \h t \ 2 G 2 )l 2 ^j . It is 
important to note that S and X are Alamouti space-time codes 
constructed at Si and S2, respectively. 

III. Average Sum-Rate of Distributed Alamouti 
Space-Time Codes with Two-Way AF Relaying 

In this section, we investigate the information-theoretic per- 
formance of the proposed DASTC with two-way AF relaying 
presented in Section[n] As mentioned above, since the channel 
model is assumed to be ergodic block fading, we can describe 
the average sum-rate as follows: 



a = E ^ l0 & det ( l2 + No (TtWf) 



h ItJ 1 !> (9) 



+ |log 2 det(j 2+ iVo(i+ 1 G , 2|fei|2 . 

where R s and R x are the covariance matrices of Gaussian 
codewords and the pre-factor 2/3 accounts for the fact that 
each source received information consisting of two symbols 



over three time slots. Due to the orthogonality of Alamouti 
codeword, it is easy to see that R s = E |SS^ j = 2V s l2 and 

{-XX t } = 2V S I 2 . Then we can rewrite © as 
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E • 



Csum — 



M 3 log 2 



I + 27 



\hi\ 2 \h 2 \ 2 
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N 2 ) 



l + 2 7 



\hi\ 2 \h 2 \ 2 



l^o| 2 (l^i 



2|/* 2 | 2 



|/» 2 | 2 +2|/ ll | 2 



(10) 



where ( [Tol l is established by the fact that det (I + AB) = 
det (J + BA). For notational simplicity, let us denote ji = 
2~/\hi\ 2 with i = 0, 1,2. Moreover, in the high SNR regime, 
i.e., 7 is sufficient large, the average sum-rate C sum given in 
([Tot , after some manipulations, can be tightly approximated 
as 
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(11) 



In order to calculate ( TPTi i. we first start our derivation by 
considering the similarity between I\ and I 2 in (fTTT >. To eval- 
uate the integral 1\ we utilize the probability density function 
(PDF) of Z derived in Lemma Q] given in the Appendix. From 
\, we have 



lnz 
l + z 



z u exv{-p lZ )K. v (az) dz (12) 



where 



a = 2 
A) = Q 



2 „ 1 2 

— , Px = — + — 

7i72 72 7i 

2 , A = ajSi 



To further simplify (fT2l . we will express lnz/(l + z) and 
exp (—az) K. v (az) in terms of the Meijer's G-function with 
the help of (TJ eq. (8.4.6.13)] and [13, eq. (6.4.23.3)] as 
follows: 
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(13) 
(14) 



where G™<"(-) is the Meijer's G-function UJ eq. (8.2.1.1)]. 
Moreover, the Meijer's G-function is a special case of the 
Fox's H-function (UJ eq. (8.3.2.21)] 
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(15) 



By combining dl3t . (fT4l . and ( Tl3T > with ( fT2l . we obtain 
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dz 



(16) 



Then, the integral 1\ can be calculated with the help of |[T4l 
eq. (2.6.2)] as follows: 
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(l + M) 
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(17) 



where -ff^'r^'.c] '^[b-d] H * s tne g enera li ze d Fox's i/-function 
|[T4l eq. (2.2.1)]. Similarly, we can obtain the closed-form 
expression for I 2 . 

We next evaluate the integral 1% by utilizing the PDF of T 
derived in Lemma [2] given in the Appendix. From d32l . we 
have 



h = J2tn f ln(l + <) (t + Q)- Z - n 
n=l g 

x 2J F 1 (3 + n,n+i;|;|±2)* (18) 

Then we exchange the variable — and as- 

sume that the relaying protocol is symmetric, i.e., 71 = 
72, yielding r\ = 0. We next express ln(l + £x) and 
(1 + x)~ a 2 Fi (a, b; c; j^jj in terms of the Meijer's G- 
function with the help of UJ] eq. (8.4.6.5)] and UJ] 
eq. (8.4.49.26)] as follows: 



ln(l + (x) = G\\ \ C 
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By combining (fT9l and (l20l with (fT8l . we have 
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The integral given in (f2Tb can be finalized in the form of Fox's 
H-function by using ESI eq. (7.811.1)] and (13 eq. (8.3.2.21)] 
yielding /3 as follows: 

(0,l),(2n,l),(l,l),(l,l) 

(2 + n,l),(n-i,l),(l,l),(0,l) 

(22) 

By substituting (fTTI i and d22l in ( fTTT i and considering the fact 
that 7 2 can be obtained in a similar form of Ii, the sum-rate 
of DASTC with two-way AF relay can be shown as 
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(23) 



IV. Numerical Results 



In this section, we provide the numerical results to verify the 
proposed two-way DASTC and the correctness of our analysis 
in two specific examples. The path loss of each link follows 
an exponential-decay model: if the distance between the two 
sources is equal to d, then ilo oc d~ a where an exponent of 
a = 4 corresponds to a typical non line-of-sight propagation 
0. Here we assume that the Si — > S2 link has unit channel 
mean power, i.e., fin. = 1- I n the first example, we consider the 
symmetric relaying protocol, i.e., ilo = fii = ^2 = 1- In the 
second example, we assume that the geometry for locations of 
three communicating terminals is co-linear where the relay is 
placed half way between the two sources, i.e., = fi 2 = 16 
and fin = 1. 

Fig. 12 and Fig. [3] illustrate the average sum-rate of DASTC 
in two-way AF relay networks versus average SNR for the 
two considered examples. As can be observed from these two 
figures, the analysis is very tight from the middle to high 
SNR regime. Specifically, from SNR=10 dB the analytical and 
simulation curves perfectly match with each other which verify 
the tightness of our approximation. 

More importantly, the proposed two-way DASTC scheme 
outperforms the conventional one-way system in terms of the 
spectral efficiency in the whole considered range of SNR. 
In particular, at SNR=30 dB, the proposed scheme enhances 
the average sum-rate to 3.8 bps/Hz in both Example 1 and 2 
compared to the conventional one. It is interesting to observe 
that with a fixed value of fin the gain is unchanged, irrespective 
of the relay's location. 

V. Conclusions 

We have proposed a DASTC scheme for two-way AF 
relay networks that circumvents the loss in spectral efficiency 
inherently occurred in conventional one-way DASTC system. 
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Average sum-rate of DASTC with co-linear relaying protocol. 



We also derive the tight approximation for the average sum- 
rate of the proposed scheme. The closed-form expression 
for approximated sum-rate is given in the form of Fox's 
H-function which readily allows us to analyze the spectral 
efficiency of the proposed scheme. The numerical results 
provided have validated our analysis. 

Appendix 

The following lemmas will be helpful in the paper. 

Lemma 1: Let 7, and jj be statistically independent and not 
necessarily identically distributed (i.n.i.d.) exponential random 
variables with hazard rate 7,; and 7,, respectively. Suppose that 
the ratio of Z takes the form 



ll + 27j 

Then, we obtain the PDF of random variable Z as 

Pz {%) — olz exp {—fiiz) [a/Co (az) + /3,/Ci (az)\ 
where a = 2./^- and ft = i- + 

V 7i7j ^' Ij 7» 



(24) 



(25) 



Proof: Let us rewrite Z in a more tractable form as 

Z=(X + Y)~ 1 (26) 

where X = 1/7,; and Y = 2/jj. Then, after some algebraic 
simplifications together with the help of [15, eq. (3.471.9)], 
the moment generating function (MGF) of X and Y can be 
expressed as, respectively 



$ x (s) = E x {cxp (-ax)} = 2 J -Id 2J- 




$y (s) =E y {exp(-st/)} 



As can be seen from (l26b . since 1/Z is the sum of two 
statistically independent random variables, by using d27l > and 
), the MGF of 1/Z can be expressed as 



V Ti7j V V 7»/ V V 7j , 



(29) 



Then, the cumulative distribution function (CDF) of 1/Z, 
Fi/ Z (x), can be shown as 



Fi/z (x) = 1-C- 1 {<S> 1/Z (s)/s}\ 1/x 



(30) 



where C 1 {.} stands for the inverse Laplace transform. From 
and (OOll together with Q6j eq. (13.2.20)], we have 



Fz (z) = 1 — az exp (— ftz) JC\ (az) (31) 
and fa = -L + Jr. By differentiating (T3Tb 



where a 
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7i7j 



with respect to z and using the fact that — ^— ■ = —K v —\ (z) — 
ZJCv 0) EG] eq. (8.486.12)], we obtain (gU which completes 
the proof. ■ 
Lemma 2: Let 70, 71, and 72 be the three i.n.i.d. expo- 
nential random variables with hazard rates 70, 71, and 72, 
respectively, then the PDF of T = 7o 7 ^^ 2 i s gi ven by 
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PT (*) (t+C) 3+ " 2-Fl ( 3 + " 
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7i72 



and 2-F1 (a, c; z) is the Gauss hypergeometric function IfTTl 
eq. (2.12.1)]. 



Proof: Since T = 7o , by defining W = 
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7172 
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Following similar steps as in Lemma Q] we obtain pw (w) as 



Pw (w) 



2w 



2/C, 



7i72 
/ 2w 
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VV7172/ V 7i72 / V V7i72, 

By substituting d34l l in (|33l and taking the derivative with 
respect to t, we obtain the PDF of T as 

2 

/ 9 / t 1 -L \ 

/ ^ \ n/n 'Yl 'Vn / 

71=1 



Pt (*) = ^ ro„ y w 2 exp 

1 

x K n 



t 1 1 

7o 7i 72 
2w 



dw 



(35) 



.a/7172, 

It is observed that the integral given in (l35T l can be simplified 
by using iTPJl eq. (6.621.3)] which results in (l32l . This 
completes the proof. ■ 
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